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Takano et al. [Phys. Rev. B 65, 140513 (2002) and unpublished] made Josephson junctions from 
single crystal whiskers of Bi2Sr2CaCu208+a crossed an angle <j>o about the c-axis. From the mesa 
structures that formed at the cross- whisker interface, they inferred a critical current density J c ((j>o)- 
Like the single crystal results of Li et al. [Phys. Rev. Lett. 83, 4160 (1999)], we show that the 
whisker data are unlikely to result from a predominantly d-wave order parameter. However, unlike 
the single crystals, these results, if correct, require the whisker c-axis transport to be coherent. 

PACS numbers: 74.72.Hs,74.50.+r,74.20.Rp 
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Recently, there have been a number of phase-sensitive 
experiments relevant to the orbital symmetry of the su- 
perconducting order parameter (OP) in the high transi- 
tion temperature Tr superconductor Bi2Sr2CaCu20s+(S 
(Bi2212).Q S S H S El II II S E3 It was claimed 
that the tricrystal experiment demonstrated a dominant 
d 3 .2_j / 2-wave OP component in Bi2212 at low tempera- 
ture T for both underdoped and overdoped samples, jlj 
However, the Pb/Bi2212 c-axis Josephson junction ex- 
periments demonstrated in many samples that at least 
a small s-wave component was present for T below 
the T c of Pb (or Nb).@, In the bicrystal c-axis 
twist experiment, |4| a dominant s-wave OP for T < T c 
was claimed.H \u\ The superb quality of the junctions 
was supported by extensive experimental and simulation 
analyses including high resolution transmission electron 
microscopy (HRTEM) studies, demonstrating that the 
junctions were atomically clean over tens of fim along 
the junction direction. The twist angle 4>o independence 
of the c-axis Josephson critical current density J c across 
the twist junction for T just below T c was interpreted in 
terms of a dominant s-wave OP for all T < T c ,0, [n| 
in apparent contradiction to the results of the tricrys- 
tal experiment. However, these experiments would be 
compatible if Bi2212 were mostly s-wave in the bulk and 
d-wave on the surface. [l2j 

Single crystal Bi2212 consists of a stack of intrinsic 
Josephson junctions, and low-T measurements of the 
critical current I c and the normal resistance R n across 
a single c-axis junction led to I c Rn values rs 1/3 the 
Ambegaokar-Baratoff (AB) result .[iO. 111 H^] This is con- 
sistent with non-metallic and incoherent c-axis transport 
in B12212, dEIIl and an s-wave OP. [hi Upon inter- 
calation with HgBr2, mesa studies revealed that I c and 
R n respectively decreased and increased by two orders of 
magnitude, but their product I c R n remained about 1/3 
of the AB value, jltj] For incoherent transport, I c for an 
s-wave or <i-wave OP is respectively proportional to the 
s-wave (1/t s ) or d-wave (I/t^) interlayer scattering rate, 
but R n oc t s for both OP's. Hence, the invariance of 
I c R n upon intercalation is strong evidence for an s-wave 



op. H3 

Recently, Takano et al. crossed two single crystal 
Bi2212 whiskers an angle 6q about the c-axis and fused 
them together, [21 |S| using a technique similar to that of 
Li et al. 4] In addition to the nominal composition of 
overdoped Bi2212, with T c rzs 80K, such whiskers usu- 
ally have a second transition at T C 2 « 105K, due to 
Bi 2 Sr2Ca2Cu3 0io +< 5 (Bi2223) contamination. 0, H For 
a 90° cross-whisker junction, a HRTEM picture re- 
vealed that this junction was uniform over 100 nm, but 
the expected periodic lattice distortion was difficult to 
discern. 9] No HRTEM pictures of other cross- whisker 
junctions were made.|ljj 

Remarkably, Takano et al. observed branch structures 
in the current-voltage (/ — V) characteristics of their 
cross- whisker junctions at 5K, indicating that the insulat- 
ing edge regions of the two whiskers had somehow fused 
into a mesa structure, Q consisting of a stack of Joseph- 
son junctions. By assuming I c of the central V = 
branch corresponded to that of the twist junction and 
that the junction area was equal to the entire whisker 
overlap, Takano et al. inferred a value for the junction 
J c for each of their 10-16 samples. U @ In sharp con- 
trast to the single crystal twist experiments, Q J c (^o) 
obtained in this way for the cross-whisker junctions var- 
ied substantially with </>o-01l3 

Nevertheless, for the 45° cross-whisker junction in a 
parallel magnetic field, a Fraunhofer-like diffraction pat- 
tern consistent with the long-junction limit was observed 
for their junction with width 36.7/im and Josephson 
length Xj ~ 3 — 4/im.@, 0] In addition, Shapiro step 
analysis of the 45° cross- whisker junction indicated that 
only first-order tunneling is present. [20j If we assume 
those results are correct, then the data for the crucially 
important region </>o ~ 45° arises from weak, first-order 
tunneling only. Since at 5K, 7 C (45°) = 0.227mA is orders 
of magnitude larger than the minimum measureablc I c , 
we use a logarithmic scale to fit the J c data. Since the 
overdoped whiskers had T c w 80 K, we take the maxi- 
mum gap value in our fits to be Ao = 22meV, consistent 
with single crystal point contact tunneling values. |2l| 



2 



Although J c (45°) ^ 0, Takano et al. nevertheless 
claimed that the strong, four-fold periodic (fro dependence 
of J c {(fro) was evidence for a "d-like" OP. 0,13 Here we 
show that such J c {4>a) behavior, if correct, is merely a 
signal of coherent tunneling in a nearly tetragonal crys- 
tal with a non-circular Fermi surface. Then, only the 
behavior of J c {(fro) f° r 4>o ~ 45°, where the Josephson 
tunneling can safely be taken to be weak and first or- 
der, is sensitive to the OP symmetry. Quantitative fits 
to the data of Takano et al. are obtained with either a 
very anisotropic s-wave gap function on the tightbind- 
ing Fermi surface, or a constant gap Ao limited to the 
extended Van Hove saddle bands, 22] but the only possi- 
bilities of a predominantly d-wave OP to fit the data are 
very unlikely. 

For weak tunneling across a cross-whisker Josephson 
junction with cross-angle <fio, J c is given by 

J c {(fro) = \4eTj2(f J (ky)F^)K$% (1) 

where we have set H = c = ks = 1, (•••) represents 
the integrations over the first Brillouin zones (BZ's) on 
each side of the junction, F w (k) = A(k, T)/[uj 2 +£ 2 (k) + 
A(k, T\ 2 ] is the anomalous Green function, A(k, T) is 
the OP with wavevector k, £(k) is the quasiparticle dis- 
persion, the u> are the Matsubara frequencies, and k rep- 
resents the wavevectors k = (k x , k y ) rotated by (fro about 
the c-axis.|I3 Here we set T c = 80 and T = 9K, suffi- 
ciently close to 5K, allowing us to include only 500 Mat- 
subara frequencies in the sum in Eq. (1). For overdoped 
cross whiskers, the normal state is rather metallic, and 
the BCS-like model, employed above, with general orbital 
OP symmetry, is a good approximation. 

For strongly incoherent tunneling, J c (<fro) / J c {0) = 
1 for all s-wave OP forms, and although for all d- 
wave forms, J c (0) is vanishingly small, J c {<fro) / Jc{0) ~ 
|cos(20q)L regardless of the details of the quasiparti- 
cle states. [ll] Neither form fits the cross- whisker data. 
0, @ However, for the quasiparticle states pictured in 
Fig. 1, a substantial fraction of coherent tunneling nec- 
essarily leads to a strong, four-fold (fr dependence to 
Jcifio) for cAo ~ 0°, 9O°,0j even for an isotropic s- 
wave OP, and J A (fro ) is mainly sensitive to the OP sym- 
metry for <fro « 45°. [ill Regardless of the tunneling 
coherence, for weak tunneling between single-component, 
nearly tetragonal d-wave superconductors, at some par- 
ticular (fro value (fro ~ 45°, the quantity inside the ab- 
solute value si gns in Eq. (1) changes sign, and hence 

Wo) = oMm 

In order to fit the whisker data, we therefore as- 
sume coherent tunneling, / J (k, k') = fo5^(k — k'), 
and take £(k) to have either the tightbinding form, 
£(k) = —t[cos(k x a) + cos(k y a)]+t'cos(k x a) cos(k„a) — u, 
where t = 306 meV, t'/t = 0.90, and fj,/t = -0.675, [H| 
or to represent the extended 'van Hove' saddle bands, 
£(k) = — 1[\ cos(k x a) — cos(k y a)\ — v\, as shown in Fig. 
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FIG. 1: Quasiparticle excitation features used in the fits. 
Solid: tightbinding Fermi surface. Dotted and dashed: Van 
Hove lines of constant energy -25 meV and -90 meV, respec- 
tively, relative to Ep. 

1. Here we take t = 500 meV and v = 2.02, so that 
the extended van Hove states never cross the Fermi en- 
ergy Ep, but lie just below it in the vicinity of the M 
points at (0,7r/a), etc. This form is a rather crude ap- 
proximation to the actual van Hove states, but serves to 
illustrate the effects of a cross-whisker junction rather 
well. Since both of these quasiparticle state forms are 
periodic, umklapp pro cesses are automatically included 
in our calculations. |23j For simplicity of notation, we set 
(fr(k) = [cos(fca;a) — cos(k y a)]N , so the ordinary d-wave 
OP is Ao</>(k), where N = 0.5 for the van Hove case, 
and N = 0.5315 for the tight-binding case, which lead 
to the maximum value of the gap Ao = 22 meV at the 
respective closest approaches to the M points in the BZ. 

A mixed d + is OP could occur below a second phase 
transition, provided that the be plane containing the pe- 
riodic lattice distortion is indeed a good mirror plane. |24| 
However, if the disorder suggested by the STM measure- 
ments on cleaved single crystal Bi2212 were present in the 
Bi2212 whiskers.flT] then there might not be any relevant 
mirror plane, and a mixed d ± s type of OP could occur 
without a second phase transition. Hence, we considered 
both OP forms. A d x 2_ y 2 + id xy OP is inconsistent with 
both the single crystal twist experiment and with the 
Pb(Nb)/Bi2212 Josephson junction experiments, 0, 0, Q 
with consequences very similar to those of the d x 2_ y 2 +is 
state. 

The d + is scenario might at first sight appear to be 
consistent with the observation of an effect claimed to be 
due to spontaneous time-reversal symmetry breaking be- 
low the pseudogap onset in underdoped Bi2212.g5J How- 
ever, in the overdoped regime of the cross-junction exper- 
iments, the angle-resolved photoemission spectroscopy 
(ARPES) measurements were only made above T c , and 
the effect was not observed, [i^ There is very strong evi- 
dence that the pseudogap exists above T c for all Bi2212 
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FIG. 2: Plot of log 10 [J c (<^ )/ < J c (90°) >] versus (f> obtained 
from Ref. [9] (solid circles). Also shown are the fits at 9K 
assuming coherent c-axis tunneling, the tight-binding £(k), 
and OP's of the ordinary d-wave, (d, long-dashed) and d-wave 
cubed, (d3, dotted) forms, along with the best d + is single 
domain fit, e = 0.f5, (solid) and the domain- averaged result 
of the best single domain d±s fit, (dot-dashed). See the text. 

dopings, and that it is distinct from the superconducting 
gap. 17, Whatever the source of the effect, 

if it were simply a property of the pseudogap, it should 
have been seen for all Bi2212 dopings. Since the effect 
was only seen in the underdoped regime, where Bi2212 is 
known to be very strongly disordered, the effect itself 
might be a combined property of the pseudogap and the 
disorder. In any event, it is unlikely to be relevant to the 
superconducting OP in the overdoped regime. 

In Fig. 2, we show our results for the best d- 
wave fits to the data, assuming the tightbinding £(k) 
form. In this and subsequent figures, we assumed co- 
herent c-axis tunneling, T = 9 K and took the value 
(J c (90°)) to be the average of the three data points 
in that vicinity. We then fit the data to the ordi- 
nary d x 2_ y 2--wave form, An0(k), a cubed <i-wave OP, 
Ao</> 3 (k), as suggested by recent ARPES measurements 
on underdoped (Pb,Bi) 2 Sr 2 CaCu208+<5, 29] a single- 
domain time-reversal symmetry broken d + is state, 
Ao[0(k) + ie}/(l + E 2 ) 1 / 2 , and a single-domain d±s state, 
Ao[0(k) ±e]/(l + |e|). For the d±s state, the best fit was 
for e = 0.25. Over the region of available data, this curve 
was nearly indistinguishable from the d + is curve with 
e = 0.15 shown in Fig. 2. In Fig. 2, we also show the re- 
sult of a multi-domain average of the d + s, d — s domains 
with e = 0.25, which includes the identical contributions 
of d ± s\d ± s domains across the cross- whisker junction, 
and the inequivalent d ± s\d =F s domain contributions. 

The best single domain d + is and d ± s fits strad- 
dled the dataless region near to 40°, and hence could be 



consistent with the data for 0.1 < e < 0.25. However, 
this time-reversal symmetry breaking d + is state has a 
vanishing J c {4>q) in the dataless regime in the vicinity of 
4>o = 40° , but not also at the crystallographically identi- 
cal cross- whisker angle 50°, for which nonvanishing data 
were available. Hence, this state would require a sin- 
gle d + is domain, in apparent contradiction with STM 
studies. 0| For a mixed d + is or d ± s state with mul- 
tiple domains, averaging over the domains always leads 
to J c (45°) = 0, as pictured in Fig. 2. We note that 
averaging over multiple d + is domains with the optimal 
e = 0.15 results in a curve that is nearly indistinguishable 
over the region of available data from the d-cubed state 
plotted in Fig. 2. 

In Fig. 3, we show the best s-wave fits using the same 
tightbinding £(k). Here we show fits to the isotropic 
s-wave OP, An, the extended-s-wave OP, An|</>(k)|, and 
several values of the highly anisotropic- s-wave OP, (An — 
A c )cj) s (k) + A c (indicated in Fig. 3 by l ext8+c-s') for 
A c = 0, 1, 1.2, 1.4, and 2 meV, as indicated. The best fits 
are for A c w 1.2 - 1.4 meV, but 1 < A c < 2 meV is ac- 
ceptable. These values are compatible with the resolution 
of recent ARPES experiments. [2^| The flat k dispersion 
away from the minimum gap position on the tightbinding 
Fermi surface is also compatible with the ARPES data. 
We note that ARPES data are complicated by the non- 
superconducting pseudoga p, which appears at T* > T c 
for all B42212 dopings. |26l l27| so that the superconduct- 
ing gap at each k of observation is merely constrained to 
be less than or equal to the gap observed by ARPES. We 
note that the curves in Figs. 2 and 3 for the ordinary 
d, ordinary s, and extended-s-wave OP's differ slightly 
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FIG. 3: Plot of log 10 [J c (<^ )/ < J c (90°) >] versus </> obtained 
from Ref. [9] (solid circles). Also shown are the fits obtained 
at 9K with the tight-binding £(k), and OP's of the isotropic 
s-wave ("s", dotted) and extended-s-wave ("ext-s", dashed) 
forms, and of the anisotropic "ext8+c-s" , (Ao— A c )<j!> 8 (k)-|-A c , 
forms with theA c values in meV of 2 (long-dashed), 1.4 (dot- 
ted), 1.2 (solid), 1, (dot-dashed), and (solid), as indicated. 
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FIG. 4: Plot of log 10 [ J c (0o )/ < J c (90°) >] versus (j) obtained 
from Ref. [9] (solid circles). Also shown are the fits obtained 
from the extended van Hove £(k), for OP's of the ordinary 
d-wave (long-dashed), ordinary s-wave (solid), extended-s- 
wave, (dotted), and the d + is forms (dot-dashed) for the 
e values indicated. For e = 1, results for the s + d state 
(short-dashed) are also shown. 

from those presented previously, since those curves 
were calculated just below T c , and these results are for 
9K<< T c . 

Finally, we present our fits relevant to the van Hove 
scenario. Here we adjusted the bandwidth t and the max- 
imum of the saddle bands — tv relative to Ep to obtain 
the best fit for the ordinary s- and <i-wave OP's. The 
values shown here, t — 500 meV and v = 2.02, are inter- 
mediate to both optimal values. In addition, we showed 
the calculations for the single domain d + is (or s + id) 
state, A [</>(k) + ie]/(l + e 2 ) 1/2 for e = 0.75,1, and 2, 
respectively. For e = 1, we also showed the results for 
the single domain s ± d state, A o [0(k) ± e]/(l + |e|). We 
note that the ordinary d-wave curve and the three curves 
with e = 0.75, 1 are inconsistent with the data, but the 
curve with e = 2 is consistent with the data. However, 
this (predominantly s-wave) d + is state exhibits a strong 
amount of time-reversal symmetry breaking, and is hence 
unlikely to be compatible with a variety of other experi- 
ments, as noted above. 

We remark that the highly anisotropic OP "ext8+c- 
s" used phenomenologically in Fig. 3 to fit the data 
could arise from a van Hove scenario, in which the domi- 
nant pairing occurs over the van Hove bands pictured in 
Fig. 1, and appears on the tight-binding Fermi surface by 
weak c oupling^ of the electronic states, as discussed else- 
where. |22j,|30j We note that a good fit could be obtained 
within the van Hove scenario using an isotropic s-wave 
OP, as shown in Fig. 4, so the physics of the generalized 
s-wave OP's in Figs. 3 and 4 need not be substantially 
different or exotic. 

Since this work was submitted for publication, it came 



to our attention that evidence exists that might cause 
one to suspect that some of the J c {4>o) values reported by 
Takano et al. might not represent the intrinsic J c {4>o) °f 
the cross- whisker junctions. For the same 90°, 75°, and 
60° cross-whisker junctions for which the mesa branch- 
ings were shown, Q the measured resistances R in the 
T region 70 - 75K« T c < T < 105K revealed significant 
complications, suggestive of the presence of more super- 
conducting Bi2223 at the 90° cross- whisker interface than 
for the 75° cross whisker, and what might be evidence for 
a non-superconducting or possibly even insulating bar- 
rier at the junction of the 60° cross whisker. Extensive 
studies of Bi2212 mesas cut from single crystals estab- 
lished that the central V — branch, with the smallest I c , 
is frequently associated with the junction closest to the 
current lead, and hence furthest from the mesa center. 
In addition, the branch with the lowest I c corresponds to 
the true critical current of the stack, [l^ These properties 
might lead to an overestimate of J c (90°), and an under- 
estimation of J c (60°), for example. Hence it is possible 
that the intrinsic 4>o dependence of J c might be weaker 
than reported. 

In summary, we have shown that the recent cross- 
whisker Josephson junction results of Takano et al., while 
different in detail from the single crystal results of Li et 
al., also render a predominantly d-wave OP form unlikely. 
As for a mixed s- and d- wave OP, there is a narrow 
window of 10-25% s-wave that would be allowed in a 
particular single domain of mixed OPs, but otherwise, 
it appears that one is forced to accept the result that 
the dominant OP is indeed s-wave, although it could be 
highly anisotropic. Although the results are compatible 
with an isotropic OP on an extended van Hove saddle 
band of quasiparticle states, if the pairing were to take 
place mainly on the tightbinding Fermi surface, the s- 
wave gap functions would have to be highly anisotropic, 
with a minimum value in the range 1-2 meV, consistent 
with ARPES experiments. |23 

In order to strengthen these conclusions, we urge that 
additional data points in the range 30° < (f>o < 60° 
be taken, and that the temperature of the measurement 
be raised up near to T c .[24| A few more junctions with 
(j)Q ss 45° are currently under study, and the prelimin ary 
results appear to be consistent with the above data.[n| 
We would also like to see other measurements to inves- 
tigate if the quasiparticle c-axis transport is indeed co- 
herent, unlike Bi2212 single crystals. We also urge that 
resistivity measurements between T c and 105K be made 
on cross- whisker junctions with 4>q ss 45°, and that alter- 
native fabrication procedures be investigated, in order to 
guarantee their uniformity. 

The author would like to thank K. Kadowaki, R. 
Kleiner, Q. Li, K. Scharnberg, J. R. Schrieffer, J. Scthna, 
and M. Tachiki for useful discussions, and especially Y. 
Takano for a preprint containing his latest data, a copy 
of Ref. [9], and extensive communications. 
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Note to be added in proof: It has ultimately come to 
our attention that an attempt to fit the whisker data of 
Takano et al. was made by Maki and Haas.JU There 
are three problems with that work. First, those authors 
assumed the Fermi surface of Bi2212 to be a cylinder 
with a circular cross-section centered about the central 
r point in the first BZ, inconsistent with ARPES ex- 
periments on Bi2212. 29| Second, they ascribed the </>o 
dependence of J c to Andreev scattering at the junction 
interfaces, which they assumed to be between super- 
conducting and normal metallic layers, [U inconsistent 
with the Shapiro step analysis and the multiple branch- 
ing behavior of the whisker current-voltage characteris- 
tics, which provide compelling evidence that the junc- 
tions are weak, first order, and of the superconducting- 
insulating variety, H, l20| as generically occurs in single 
crystal Bi2212.[lH Q Third, their formulas (3) and (6) 
are inconsistent with those obtained for their model from 
the exact expression, Eq. (18) of Ref. [23, obtained for 
the c-axis critical current for coherent tunneling between 
two layered superconductors. Hence, the J c {4>o) obtained 
by Maki and Haas is both mathematically imprecise and 
physically very unlikely to apply to Bi2212 whiskers. 
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